Remarks on (1 — q) expansion and factorization approximation in the Tsallis 

nonextensive statistical mechanics 
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The validity of (1-q) expansion and factorization approximations are analysed in the framework 
of Tsallis statistics. We employ exact expressions for classical independent systems (harmonic 
oscillators) by considering the unnormalized and normalized constrainsts. We show that these 
approxiamtions can not be accurate in the analysis of systems with many degrees of freedom. 
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I-Introduction 
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Ever since the presentation by Tsallis [jl], |], |) of a new 
possible generalization of the statistical mechanics (Tsal- 
lis statistics) based on the nonextensive entropy, 



Sq — k- 



Trp Q 
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(1) 



a large number of investigations has been developed con- 
cerning this subject B. These investigations are basi- 
cally employed in the discussion of aspects related to 
nonextensive phenomena, such as, Levy-type anomalous 



self-gravitating 
||, [l(J , cosmic 
peculiar vcloci- 



superdiffusion |5|, Eulcr turbulence) £], 

systems and related the mes f | pi 3, 

background radiation (T| g [u|, 
ties in galaxies clctron-phonon interaction |L6l, and 
ferrofluid-like systems ]l7|]. In addition, some impor- 
tant methods of the usual Boltzmann-Gibbs statistics 
have been generalized in order to incorporate the Tsal- 
lis framework, e. linear response theory fl8|| , Green 
function theory |l9) and path integral^]. However, ex- 
act calculations are generally difficult to be performed 
in the context of Tsallis statistics. Motivated by these 
difficulties some approximated methods, such as (1 — q) 
expansion | H] , |i"4| , pl| , p2[ , factorization approximation 
p3[ , perturbative expansion p4[ , generalized Bogoliubov 
inequality (2l|, ^5|, and semi-classical expansion p^, have 
been developed. The two first methods mentioned above 
have been employed largely in the analysis of black body 
14[ and in other many independent par- 
but without a careful 



radiation |TT], 13 

tide systems |29J, gCj, |3J, g2| 

analysis of the validity of these methods. Thus, a detailed 
discussion of the (1 — q) expansion and the factorization 
approximation plays a special role in this scenario. In 
this direction, Ref. J33[ contains a discussion compar- 



ing these approximations in the context of the quantum 
gases. However, some important aspects of the compari- 
son are still lacking, in particular the importance of the 
degrees of freedom N. This work is addressed to analyze 
how precise the (1—q) expansion and the factorization ap- 
proximation are by considering N arbitrary. To perform 
this study it is convenient to consider a solvable model 
in order to carefully understand the degree of accuracy 
of these approximations. Furthermore, the chosen model 
must contain important features of other more realistic 
ones. This is just the case of a set of harmonic oscilla- 
tors. In this direction, to analyze the (1 — q) expansion 
and the factorization approximation, we focus our dis- 
cussion on the classical Tsallis statistics by considering 
N one dimesional harmonic oscillators, i. e., we employ 
the Hamiltonian 



H = 



N 

E 

n=l 



2m 



(2) 



where tt^ + imw^x 2 , is the Hamiltonian to the rt-th har- 
monic oscillator, being m the particle mass and u n the 
harmonic oscillator frequency. 

This work is organized as follows. In Sec. II, we inves- 
tigate both approximations by using unnormalized con- 
straints. The analysis of these approximations employing 
normalized constraints is done in Sec. Ill and the con- 
clusions are presented in Sec. IV. 
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FIG. 1: In order to illustrate the differences among Eqs. (|l2[), 
(^) and ( |l6| ) we ploted them versus N for q — 0.99 and 
q = 0.9999. 



II-Exact versus approximate calculations and 
unnormalized constraints 



In the Tsallis statistics with unnormalized constraints, 
the nonextensive canonical distribution is given by 



Pq = ^-[l-(l-q)fiH] 1/{1 - q) , 



(3) 



where 



N 

z z = J II d Pndx n [l - (1 - -z)/3ff] 1/(1 - 



9) 



(4) 



is the partition function with q 6 R and /3 = 1/T 
(with fc = 1). Eq.(0) is obtained maximizing the Tsallis 



entropy]]], ||, Eq. ([!]), subject to the constraints 

C/, = J I] d Pn dx n p« i? (5) 



and 



/ ]^[ dp„dx n p 9 = 1 , 



(6) 



with C/ g being the generalized internal energy. In this 
context, it is employed 



r N 

A q = (A) q = / II d Pn dx n P q A 
J n=l 



(7) 



as the generalized mean value of the classical function 
A(jp, x) and thermodynamical functions like free energy 
and specific heat are defined as F q — (Z^ 1 — 1)/[(1 — 
q)(3] , and C q = dU q /dT. Furthermore, the Legendre 
structure of the Tsallis statistics is preserved ||, |3J, ^5) ; 
thus, F q = U q -T S q , S q = -dF q /dT, etc. In Eq. 
d) we assumed that 1 - (1 - q) (3H > 0. When this 
condition is not satisfied we have a cut-off, i.e., p(H) = 
if 1 — (1 — q) j3H < 0. Thus, when a classical partition 
function is calculated, the integration limits on the phase 
space are obtained from the condition 1 — (1 — q) (3H > 0. 

The partition function (Q) for the N harmonic oscilla- 
tor system is 



Z q = Jf[ dp n dx n [1 - (1 - q)/3 



N 



vfc=l 



1/(1-9) 



(8) 



To calculate this integral we introduce the variables: 
Vn = [(1 - q)muj 2 [3/2] 1 / 2 x n and y N+n = [(1 - 
q)(3/{2m)} 1 / 2 p n , where n = 1,2,3,...,N. In terms of 



these variables, Z q becomes 



N 



n 

v n— 1 



(1 - q)ui n P 



r 2N / IN \ 

[ n dynil-a^yt) 

J n=l \ fe=l / 



1/(1-9) 



(9) 



with a — sign(l — q). For sufficiently large N this integral 
diverges when a = —1, i.e. when q > 1 + l/N. This fact 
indicates that the Tsallis statistics with q > 1 can not 
be employed for large number of subsystems. At this 
moment, we restrict our discussion to the case q < 1. By 
using hyperspherical coordinates with u = (X}«=i Vn) 1 ^ 2 
and performing the integral over the angular variables, 



3 



we obtain 



Z„ 




n 



2N 



(1 - q)uj n P_ 

i JV-1 / 1 \l/(l-g) 
Q11H (1 — u) ' v 



(10) 



By substituting the expression for the solid angle) 36 
— 2ir N /T(N), and employing integral representa- 
tion of Euler beta function |3q], we verify that 



' N 

n 



2tt 



(1 - q)u n (3. 



' 1-9 



1 



2-g 



-1 !f r 



n 



2tt 



+ l + iV 



(1 - q)w n (3 



(11) 



where (a) n = Jlfe=o ( a + ^) * s ^ ne Pochhammer symbol. 

In the following analysis, one can consider to n = lu 
without loss of generality. In fact, if someone is interested 
in recovering the case w„ ^ to, it is sufficient to replace 



with n^=i w n m t ne following equations. 



To investigate the (1 — q) expansion and factoriza- 
tion approximation we can employ any thermodynamical 
function. In this work, we use the internal energy. As 
exposed below Eq. (Q), the generalized internal energy 
can be obtained from the identity U q — F q — TdF q /dT, 
so it is given by 



U a 




2tt 



q)ujf3 



N r 



1-3 



-i -l 



JV 



1-9 



.(12) 



In order to analyze the (1 — q) expansion, we consider a 
series of Eq. ( [l2j ) in powers of (1 — q), 



N 


"2tT 


{l-q)N ■ 


1 







\--{l-q) 2 N{N+l) 



12 



(l-q) 3 N(N + 1){2N + 1) + 



(13) 



where the upper index factor refers to the quantities ob- 
tained with the factorization approach and Z q repre- 
sents the partition function of the subsystem labeled by 
i, and N is the number of subsystems. In the present 
study, each Z q is the partition function of one harmonic 
oscillator. Thus, the partition function Z^ actor can be 
obtained from (TnT) by using N = 1, i. e. 



r? factor 



N 

n 



2tt 



\ (2 - q)up 



2tt 



(2 - q)u>P 



N 



(15) 



where we employed = 2tt/((2 — q)u>(3) and w„ = lu. 
The internal energy in this case is given by 



jjfactor 



2tt 



(2 - q)u0 



(l-q)N 



(16) 



Afte r p erforming an expansion in powers of (1 — q), 
Eq.(|l5|) can be written as 



tj factor 
u q 



N 
1 
N 
J 



2tt 



(2 - q)u>0 



2vr 



u/3 



(l-q)N 



l-(l- q ) 2 N--(l-q) 3 N + 



.(17) 



The comparison between ( 13 ) and ( [T7j ) shows that (1 — q) 
and N dependences in U q and U^ ac come out to be dif- 
ferent. Similar differences are also present in the other 
thermodynamical functions. Furthermore, it is impor- 
tant to remark that the N dependence on Eqs. ( |l2| ) and 
( |l6| ) differs drastically for large N (see Fig.(l)). In addi- 
tion, we also see in Fig.(l) that Eq.jT^) and Eq.([l3]) are 
closed for the values of q and N when (1 — q)N <C 1 is 
verified. However, when this condition is not satisfied, 
these equations lead to a significative deviation from the 
exact result, Eq. ( |l2| ) (see Fig.(l)). Furthermore, the 
error due to the approximation is more pronounceable 
for the factorization approximation than for the (1 — q) 
expansion. Thus, U q factor is not also an accurate approx- 
imation for U q when we consider an arbitrary q and a 
very large N. 



Since the nth term of this series essentially contains the 
factor [(1 — q)N] n , the convergence of the series is im- 
proved if (1 — q)N <C 1. In other words, the above ex- 
pansion seems to be reasonable for arbitrary f3, m and u> 
only if (1 - q)N < 1. 

Let us now analyze the possibility of factorization of 
the partition function of a system with N subsystems. In 
other words, we want to understand when the partition 
function can be well approximated by 



N 



z factor = -Q z (i 



(14) 



III- Exact versus approximate calculations and 
normalized constraints 

In this section, we analyze the differences among the 
factorization approximation, (1— q) expansion and the ex- 
act calculation in the normalized approach of the Tsallis 
statistics. The nonextensive canonical distribution p" 
which emerges from the maximum entropy principle [J 
and subject to the normalized constraints^ 



jjnor 

9 — 



Tr (p- 



H 



i=l 



Tr (pj") 



and 



Trp- 



= 1 



(18) 
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N 



for the exact case. In particular, Eq. ( pl[ ) expressed in 
terms of the effective temperature T e // = Tr(p™ or ) 9 //3 
(for a discussion about effective temperature and La- 
grange parameters see [|[ |t], [||) is U q , nor = NT e ff. 

Let us now evaluate the internal energy by considering 
(1 — q) expansion and the factorization approximations 
in order to perform the analysis of these approximations 
in the normalized approach. We start by expanding Eq. 
( pl| ) in series of (1 — q) as follows: 



u =>L 

^q.nor — n 



2tt 

u/3 



1-(1-<J)N l+(l-g)JV 

[1 + (1 - q)N] i-d-9)« 



l--(l-q) 2 N(N + l) 



- —(l-q) 3 N(N + l)(m-l) + 



(22) 



The internal energy in the factorization approach can be 
obtained employing similar calculation as the previous 
one in Sec. II and it is given by 



jt factor 
q,nor 



N 
J 



2vr 



(2 - q)u/3 



l + (l-g)JV 

X [1 + (1 - q)N] 



N 

J 



2tt 

w/3 



l-(l-q)JV l + (l-g)W 

[1 + (1 - q)N] r= ^=^ r 



1 - (1 - q) 2 N + -(1 - q) 3 iV(2iV - 1) + 



(23) 



FIG. 2: In this figure, we show the differences among Eqs. 
(|l|), @ and (||) versus TV for g = 0.99 and q = 0.9999. 



is given by 

nor 1 

/9„ = — 

' 4 ynor 



l-(l-q)-^(H-U, 



Trp 



nor\ 
q I 



1/(1-9) 



where Z™ or is the partition function 



Z q = Tr 



l-(l-9)^|(^-^ n ° r ) 



1 1/(1-9) 



(19) 



(20) 



In general, the calculations in the normalized Tsallis 
statistics may be obtained from the unnormalized one||. 
In this way, employing Eqs. (|l8|) and (19) with the 



Hamiltonian given by (|2j) and taking into account Eq. 
(O), we obtain that 



U, 



q,nor 



N 
J 



2tt 



1 N 



(1 - q)ujp 



n-ll l-(l-q)JV 



l + (l-g)JV 

x [1 + (1 - q)AT] i-(i-9>jv 



(21) 



The dependence on N and (1 — q) in the expansions (gS 
and (^3|) is notoriously different. In fact, it essentially re- 
members the previous one performed in Sec. II. Thus, the 
conclusions about how useful (accurate) the (1— (?) expan- 
sion and the factorization approximations are when de- 
scribing the behavior of the thermodynamical functions 
in normalized cases are similar to the ones presented in 
Sec. II, i. e. , the approximations seem reasonable only 
for (1 - q)N < 1. The N dependence on Eqs. @, @ 
and ( p3| ) is also illustrated in Fig. (2). Moreover, as in 
the case of the unnormalized Tsallis statistics, the (1 — q) 
expansion gives a better agremment with the exact result 
than the factorization approximation. 

It is important to remark that similar results to those 
obtained in Sees. II and III are also valid when q > 1. 
In this case, a further care is necessary as the inequality 
q < 1 + 1/N must be obeyed so the partition function (Q) 
remains finite. 



IV-Conclusion 

The previous discussions demonstrated that the (1 — q) 
expansion and the factorization approximation are not 
useful for a system with N harmonic oscillators when 



5 



someone employs an arbitrary q and a very large N. In 
particular, for systems where the number of oscillator is 
taken as the Avogradro number, i.e., N w 10 23 , q ~ 1 is 
implied, so it is consistent with the discussion presented 
above. Similar conclusions can be directly verified for 
the classical ideal gas, since the calculation details re- 
lated to it are similar to that performed here for a set 
harmonic oscillators. We would like to stress that the 
system discussed here is usually employed as a first ap- 
proximation in many contexts and consequently it has 
received a special attention in the studies based on the 
nonextensive Tsallis statistics. In conclusion, the results 
presented here strongly suggest that (1 — q) expansion 



and factorization approximation are not good approxi- 
mating methods when q is generic and N is very large; 
thus, we should be very careful when employing these 
approximations. 
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